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In 2007, a very bright radio pulse was identified in the archival data of the Parkes Telescope in Australia, marking the beginning 
of a new research branch in astrophysics. In 2013, this kind of millisecond bursts with extremely high brightness temperature 
takes a unified name, fast radio burst (FRB). Over the first few years, FRBs seemed very mysterious because the sample of known 
events was limited. With the improvement of instruments over the last five years, hundreds of new FRBs have been discovered. 
The field is now undergoing a revolution and understanding of FRB has rapidly increased as new observational data increasingly 
accumulates. In this review, we will summarize the basic physics of FRBs and discuss the current research progress in this area. 
We have tried to cover a wide range of FRB topics, including the observational property, propagation effect, population study, 
radiation mechanism, source model, and application in cosmology. A framework based on the latest observational facts is now 
under construction. In the near future, this exciting field is expected to make significant breakthroughs. 
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1 Introduction 


Fast Radio Bursts (FRBs) are transient radio pulses of mil- 
lisecond durations that flash randomly in the sky. They were 
first discovered by the Parkes Telescope in 2007 [1]. De- 
spite being confused with rotating radio transients [2] and 
“perytons” [3] in the beginning, FRBs were finally deemed 
to be new astrophysical signals, thanks to a new sample of 
FRBs identified in 2013 [4]. Since then, the scientific com- 
munity has been increasingly focusing on them and theorist- 
s have been investigating the underlying physics. Knowl- 
edge regarding FRBs has taken a great leap forward since 
the landmark discovery of the first repeating FRB 121102 
[5,6] and further the identification of its host galaxy [7-9]. 
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In the following years, an increasing number of events 
have been collected [10], benefiting from the improved ob- 
servational techniques and new facilities. At the time of 
writing about twenty repeaters have been reported mainly 
by the Canadian Hydrogen Intensity Mapping Experimen- 
t (CHIME) [11-13] which has a large field of view and is 
extremely powerful in finding new FRBs. 
metric arrays, represented by Australian Square Kilometre 
Array Pathfinder (ASKAP), can accurately localize individ- 
ual bursts [14, 15]. The recent two breakthroughs in 2020 
are even more exciting: the periodical activity has first been 
found for the well-known repeater FRB 180916 [16] and later 
for FRB 121102 [17, 18]. Moreover, the association of FRB 
200428 with short X-ray bursts (XRBs) from a Galactic mag- 
netar SGR 1935-2154 has been confirmed [19-25], shedding 
new light on the origin of FRBs. 


The interfero- 
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As we deduce from the above milestones, the field of FRB 
is currently flourishing. The readers may refer to some previ- 
ous good reviews [26-29] for the progress till 2019. As this 
field is evolving fast and our understanding is being constant- 
ly updated, we feel it is a proper time for the latest review to 
appear. At the time of writing, we noticed a newly-published, 
well-written review article on the FRB radiation mechanism 
[30], and a brief introduction of recent observational progress 
[31]. With better data resolution and a larger sample of FRBs, 
we are now reaching a comprehensive understanding of this 
phenomenon. First, we get to know some individual bursts 
better, including their pulse structure, spectrum, polarization 
properties, and radiation mechanism. Second, the population 
research is also possible and we could study the classification, 
event rate and a few statistical properties of FRBs. Last, there 
is a growing trend that FRBs are treated as powerful probes of 
fundamental physics and cosmology. In this review, we will 
discuss all these aspects on the basis of the latest discoveries. 

This review is organized as follows. In Section 2 we in- 
troduce the basic observational facts of FRBs. Section 3 & 4 
present the recent progress on statistical properties and pop- 
ulation study of FRBs. The possible radiation mechanisms 
and source models of FRBs are discussed in Section 5 & 6 
respectively. The application of FRBs in cosmology is ex- 
plored in Section 7. Finally, we provide a summary and some 
thoughts regarding future prospects in Section 8. 


2 General Observational Properties 


There is no strict definition for FRBs right now. Generally, 
an FRB is defined as a radio pulse with duration T ~ ms and 
extremely high brightness temperature Tg, which can be ob- 
tained by equaling the observed intensity 7, with blackbody 
luminosity at the Rayleigh-Jeans end, 


I, = 2kTg/ A2, (1) 


where k is Boltzmann constant and A = c/v is the radio wave- 
length. Z, can be related to flux density F, as I, ~ F,/Q, 
where the observed solid angle is determined by the source 
size and angular distance, i.e., Q ~ n(cTy /d* [30]. Thus, 
the brightness temperature of an FRB with typical scaling is 


Tg ~ Fyd [2nk(vTY. 


= ax (Sz) C.) (is) 0 


which is several orders of magnitude higher than that of the 
normal pulsar radio emission. 


Radio waves can be easily affected by the intervening 
medium between the source and Earth observer. First, the 
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radio pulse is dispersed in the plasma, i.e., the high frequen- 
cy photons arrive earlier than the low energy ones. Second, 
the pulse can be broadened due to scattering with an extend- 
ed screen. Third, diffraction and refraction by turbulent gas 
make the radio signal scintillate occasionally. Fourth, these 
refractive plasma could lead to lensing effect. Fifth, some- 
times the absorption of radio emission is important if the plas- 
ma is dense enough. Last but not the least, the polarization 
property can be changed as radio wave propagates in the mag- 
netized plasma. We will introduce these propagation effects 
of FRBs respectively in the following subsections. 


2.1 Dispersion 
The dispersion relation in a cold plasma is 
w = w, +k°c?, where wp = N4nn,e? [m,, (3) 


so the group velocity of the electromagnetic (EM) wave is 


Vg < = c4J1 - (wp/w)?. (4) 


The delay time of a signal traveling through the plasma is 


Ja 
Il 


- [60 8G) + a(S) +-}o 


In radio astronomy the dispersion measure (DM) is defined 
as the electron number density integrated along the traveled 
path DM = f n,dl/(1 + z) where z is the redshift, so that the 
delayed arrival time of two photons with frequencies v1, v2 
(vı < v2) can be expressed as [28,32], 
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Figure 1 shows the dynamic spectrum (“waterfall” plot) 
of the first-discovered event, FRB 010724 (also know as 
“Lorimer burst") [1]. The dispersive “sweep” clearly indi- 
cates that high energy photons arrive earlier than low energy 
ones. This burst has a total DM of about 375 pc cm^?, which 
exceeds the DM budget that could be contributed by electrons 
in the Milky Way, suggesting an extragalactic origin. 
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Figure 1 Dynamic spectrum of the Lorimer burst. (Adapted from [1]. 
Reprinted with permission from AAAS.) 


Generally, the total DM of an FRB at redshift z can be sep- 
arated as 


DM OS F DMsource 
DM = DMyw + DMicm + —— x: zo. d 
Z 


where the subscripts MW, IGM, host, source refer to the con- 
tributions by the plasma of Milky Way, intergalactic medi- 
um, host galaxy and source local environment respectively. 
DMywvw can be obtained from Galactic electron density mod- 
els such as NE2001 [33] and YMW16 [34]. It is also useful to 
define an excess of DM with respect to the Milky Way value, 
i.e., 

DMhost + DMource 


DMg =DM- DMww = DMicm + 
1+z 


(8) 


The average (DMjgm) as a function of redshift is given by 


[32]: 
_ 3cO Ho < F(z’) " 
Mow = or [mede ©) 


where Ho is the Hubble constant, Oy is the cosmic baryon 
mass density fraction, E(z) = HA(z)/Ho and F(z) = (1 + 
z)fiaM (2) fe (z) with fiom being the fraction of the baryon mass 
in the IGM. Further, f. = YyXe.u(z) + 3 YueX«ne (2). where 
Yu, Yue are the mass fractions of hydrogen and helium, and 
XeH, Xe ge are the ionization fractions of intergalactic hydro- 
gen and helium, respectively. It is worth mentioning that the 
DMyjyay term has been applied widely in cosmology, which 
will be discussed in detail later in Section 7. 


2.2 Pulse Width and Scattering Effect 


The pulse width is commonly defined as the full width at half 
maximum. The observed FRB width is a combined result of 
the intrinsic width W;, instrumental broadening and, propa- 
gation effect [28], 


= 2 2 2 2 
W= a WF + fzmp + Alpy + Mpeg 


tT, (10) 
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where fsamp is the data sampling interval and r is the scatter- 
ing timescale. The dispersion smearing due to finite channel 
bandwidth Av of the receiver is [35] 


Atpm = 8.3 us DMAvunzYai- cm) 


Also, an error of trial DM (i.e., DM.) with respect to the true 
value in dedispersion could cause additional smearing, which 
is expressed as Atpmerr = Atpw(DMer,/DM) [35]. 

A radio photon can be easily scattered by the particles on 
its path. As aresult, the moving direction of this photon could 
be changed to some extent. For an FRB event, a portion of its 
photons could be scattered and travel a longer way then arrive 
later than the unscattered ones, giving rise to a tail feature of 
the pulse profile. Usually, this tail appears as an exponen- 
tial decay, which can be well modeled by the scattering with 
a thin extended screen. The scatter-broadening time scales 
with frequency as T œ v * and œ ~ 4.0 for the above screen 
[36]. However, if a Kolmogorov spectrum of scattering medi- 
um is assumed, the index could be a ~ 4.4 if the minimum 
turbulence scale is smaller than the diffractive length scale 
[35, 37]. 

Practically, the scattering timescale of an FRB is deter- 
mined by fitting its pulse profile. First we need to make a few 
assumptions on its intrinsic pulse shape (usually assumed a 
Gaussian shape), the receiver noise and instrumental smear- 
ing. Distinctions in assumptions lead to a diversity in FRB 
signal models. Then these models are applied to each indi- 
vidual FRB through a Markov-Chain Monte-Carlo (MCMC) 
fitting of the pulse flux. Last, the goodness of fitting and the 
preferred model is determined by the Bayes Information Cri- 
terion. Ravi (2019) collected a sample of 17 Parkes FRBs and 
studied their scattering properties [38]. He found the typical 
scattering timescale r ranges from submilliseconds to tens of 
milliseconds. Further, FRBs are underscattered in the plane 
of t — DMg, compared with the 7 — DM relation of Galac- 
tic pulsars [39-41]. These have been confirmed by a latest 
analysis on a sample of 33 ASKAP FRBs [42]. 


2.3 Scintillation 


The scintillation of FRBs is akin to the twinkling of stars. 
However, the latter is caused by turbulent atmosphere of the 
Earth, while the former is due to turbulent plasma medium. 
Once the turbulent intervening plasma moves perpendicular- 
ly to our line of sight (LOS) with a velocity V, the interfer- 
ence fringes and diffraction fringes would vary according- 
ly and observer on Earth could be in the position of either 
maximum or minimum intensity. This leads to a variation of 
observed flux density, depending strongly on the frequency. 
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Also, the variability of scintillation pattern depends on turbu- 
lence spectrum and the relative motion between the medium 
and observer. 

The physics of pulsar scintillation has been well reviewed 
before [43, 44]. Random phase fluctuations $(x, y) are gen- 
erated when the plane waves pass through a turbulent scat- 
tering screen. The diffractive length scale rai is defined as 
the transverse radius at which the root mean square phase d- 
ifference is 1 rad [44]. If raim is larger than the Fresnel scale 
rg, Which is expressed as rp = VAD/2z with D being the 
distance between the screen and observer, the perturbations 
to wavefronts are weak and this is called the weak scattering 
regime. The observed scintillation timescale is the Fresnel 
timescale tsin ~ tp = rp/V. In the opposite strong scat- 
tering regime (ra < rr), diffractive interstellar scintillation 
(DISS) plays an important role and the scintillation timescale 
iS fscint = rgig/V. Detailed discussion of the dependence of 
DISS timescale on velocities of the source, scattering medi- 
um and observer was given by Cordes & Rickett (1998) [45]. 
They found that the scintillation bandwidth Aygcint is relat- 
ed to scattering T by 27AVecintT = C, where C, is a factor 
of order unity. Combined with the previous discussions on 
T, we get AVscint © Y^. Scintillation phenomena that have 
been observed in a few cases like FRB 150807 and FRB 
121102 match the predictions well [46, 47]. Note that the 
spectral structure of the recent Galactic FRB 200428 could 
be potentially explained by scintillation [48]. Besides, the re- 
fractive interstellar scintillation (RISS) may occur on a larg- 
er scale Fef = "ras, of which the scintillation timescale 
tscint ~ Fret /V could be much longer than the FRB burst dura- 
tion [6]. 


2.4 Plasma Lensing 


Ionized plasma is refractive that radio waves are deflected 
while passing through. The plasma refractive index is given 
by np = ./1 - w/w. In many ways plasma lensing is anal- 
ogous to gravitational lensing, however, the differences be- 
tween them are also obvious. In general, gravitational lenses 
are converging and achromatic [49], while plasma lenses are 
diverging and highly chromatic [50]. If the diverging plasma 
lens is on the axis between the source and the observer, the 
observed flux will be in its minimum. However, this lens can 
also cause large magnification (caustic spots) if it lies offset. 
Therefore the variability of light curves is enriched due to the 
source-lens-observer geometry. 

The mathematical prescription of plasma lensing is very 
similar to that of gravitational lensing. The major difference 
lies in the formulae of effective deflection potential and a de- 
tailed comparison can be found in Wagner & Er (2020) [51]. 
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The electron density distribution on the lens plane will in- 
fluence the deflection angle to a great extent and differen- 
t lens models have been developed. The most commonly- 
adopted model is one-dimension Gaussian lens [50, 52] and 
more complicated light curves can be realized if other lens 
models are assumed [53, 54]. 

Plasma lensing effect could be relevant for FRBs in various 
aspects. It can introduce an additional frequency-dependent 
delay of the arrival time so that the inferred DM should be ex- 
amined carefully [55]. The impact of plasma lensing on the 
observed FRB luminosity and spectrum has been discussed 
[52]. Moreover, the complex time-frequency pulse structure 
of FRB 121102 could be possibly explained by plasma lens- 
ing [56]. The downward frequency drift in a sequence of 
sub-bursts (*sad trombone" effect) can be attributed to plas- 
ma lensing. However, upward drifting is also expected in this 
scenario but has not been observed in the majority of FRB- 
s. Even though the sign of upward drifting appears for FR- 
B 190611 [57] and FRB 200428 [20], it happens much less 
frequently than downward drifting. Note that several inves- 
tigations have been done to explain the sad trombone effect 
in distinct scenarios. Within the framework of coherent cur- 
vature emission in the magnetosphere, Wang et al. (2019) 
proposed a generic geometrical model that EM waves of dif- 
ferent frequencies are emitted by bunched particles streaming 
along field lines of different curvature radii [58]. Alternative- 
ly, Metzger et al. (2019) suggested that the sad trombone ef- 
fect could be attributed to the deceleration of the blast wave, 
assuming the FRB mechanism is synchrotron maser emission 
[59]. A detailed discussion of this model can be found in Sec- 
tion 6.1. 


2.5 Polarization and Faraday Rotation Effect 


Polarization refers to the behaviour of the electric field vec- 
tor E with time. Owing to the variety of the trajectory that 
E draws on the plane perpendicular to wave propagation, we 
have linearly, circularly, elliptically polarized or totally un- 
polarized light. Suppose an EM wave propagates in the same 
z-axis direction with the magnetic field of magnetized plas- 
ma (B = Boe-), then E oscillates in the x — y plane and the 
equation of motion for an electron is 


Mep = —eE — ey X (Boe.). (12) 


For the plane wave E = Ege *'-^^, and the electric polariza- 


tion intensity P = —n,er, = (e— €))E, we have rp = -E9 
Substituting into Eq. (12), we get 
2 
Mew . €— & 
| (e — ey) * e| E, + iwBo 2 E,-0 
Nee ES 
Mew . €— & 
(e — &) * e| E, - io Bo E, =0. (13) 
Nee : z 
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The determinant of coefficient should be zero, therefore 


Mew E e- ey 
(e—- e) +e] =| wByo—— (14) 
E ne 
Since the refractive index n? = €/&, we obtain 
1 oe E (15) 
1n w + oO, i 


where +, — represents left and right handed wave respective- 
ly and Q, = eBo/mec is the cyclotron frequency. A linearly 
polarized wave can be decomposed into left- and right-hand 
circularly polarized components 


E = E,(e, — ie)e ^? + E (e, tiei, (16) 


The Faraday rotation effect is the rotation of polarization 
direction of linearly polarized light under the influence of a 
magnetic field. Suppose that the magnetized plasma extends 
from z = 0 to L and without loss of generality we assume the 
EM wave is initially linearly polarized in the x direction, i.e., 
E(z = 0) = Ege "e, so E, = E. = (1/2)Ep. The electric 
field vector at the exit of the medium is 


E(z = L) 


1 . i : 
5 Eo [e (e, — ie) + e" (e, + ie) e" 


1 M 
- KL, EL CL JL 
5 Boe Vr IG t t ee, rie — e e; | 


ky — k- k, — k- TS 
= Eo [cos i Le, ^ sin — Le, e Loon 
(17) 
Now this wave is linearly polarized at angle 
k, — k. w L 
Y= x L- c0. -n25 (18) 
In the limit w >> Q,, Eq. (15) can be approximated as 
2 
&l-.————, 19 
"m 2 0? + wQ, (19) 


thus in general case, the polarization direction after traveling 
an infinitesimal length d/ will rotate: 


2 
wdl cx Odi 2ne3 
d¥ = oe —n)s cu = nee neBodl. (20) 
Therefore, the observed polarization angle (PA) 
ez L 
Pops(A) = Yo + AV = Yo + ——— f ne(DB (Ddl, (21) 
2nm2c^ Jo 


where Yo is the initial PA and Bj is the component of the 
magnetic field along the LOS. In terms of the rotation mea- 
sure (RM), 


3 


L 
RM if neDB\Ddl, (22) 
0 


2am2ct* 
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Vos, can be written as 


WP ops(A) = Yo + RM. (23) 


For cosmological sources, RM value is usually scaled as 


L -3 
pM] =| = ost2 | rele "BIG i 
0 


m? (14 2? d 


A positive RM means the magnetic field direction is towards 
the observer and inversely it can be negative. 

Orthogonal feeds of radio telescopes can measure the ful- 
] Stokes parameters (/, Q, U, V) and PA can be determined 
by V, = 0.5tan^! (U/Q). Then one can fit the relation with 
wavelength in Eq. (23) and estimate RM. More precise RM 
value can be obtained by means of RMFIT, RM Synthesis or 
QU-fitting [60-62]. Further, if DM and RM originate from 
the same plasma medium, we can estimate its magnetic field 
strength as [63] 


RM DM \' 
~ 1232 : —| 4G.Q5) 
rad m~? JA pc cm? 


The polarization properties of FRBs vary from case to 
case. The repeaters like FRB 121102 and FRB 180916 are 
nearly 100% linearly polarized and have a flat PA curve (i.e., 
PA does not vary obviously with time) [12, 47, 64-66]. Note 
that the RM of FRB 121102 is the highest among all FRB- 
s (^ 10? rad/m?) and is decreasing with time [64, 67]. The 
non-repeaters like FRB 181112 and FRB 190102 are usual- 
ly partial linearly or circularly polarized and their PAs vary 
significantly with time (named “PA swing") [57, 68]. The 
polarization fractions are found to evolve with time for non- 
repeaters [57]. However, it remains questionable whether the 
difference in polarization favors the classification of FRBs 
by repeating or not. In particular, a recent observation chal- 
lenges this viewpoint, i.e., the repeating FRB 180301 shows 
diverse PA swings [69]. It is possible that with the enlarging 
of the FRB sample in the future, repeaters and non-repeaters 
are not so different in polarization. 


RM 
(B) = Sapa 
0.812DM 


2.6 Absorption Processes 


The most relevant absorption process in radio band is free- 
free absorption. Free electrons jump to higher energy states 
after absorbing photons and the optical depth of this process 
depends on the cross section and electron density. General- 
ly the IGM is sparse and the free-free absorption in IGM is 
unimportant. Still, significant absorption could possibly hap- 
pen in the vicinity of the FRB source, where the plasma is 
dense enough. For instance, the FRB progenitor may locate 
in pulsar wind nebulae where there are numerous electrons 
[59, 70-72]. In this case, synchrotron self-absorption (SSA) 
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may also play an important role [73, 74]. FRBs can escape 
only if the optical depths tg and Tssa are less than unity. 

Early in 2020 the lowest frequency of FRB detection- 
s is still near 330 MHz for FRB 180916 [75, 76] and FR- 
B 200125A [77]. Very recently, searches conducted at even 
lower frequencies (i.e., 150 MHz by LOFAR) found a few 
new bursts for FRB 180916 [78, 79]. This means that the 
low-energy spectrum cutoff lies below ~ 100 MHz. It is un- 
clear whether the low-energy cutoff is caused by absorption 
or intrinsic FRB radiation spectrum. However, a general con- 
straint on the intrinsic spectrum can be given for particular 
bursts [75, 80]. 


2.7 Host galaxies 


The host galaxies of FRBs are crucial for understanding the 
environment and origin of FRBs. The first localized even- 
t is repeating FRB 121102 [7]. 
metallicity, dwarf galaxy with stellar mass (1.4 + 0.7) x 105 
Mo. A radio continuum spectrum has also been found in this 
galaxy. The second localized repeating FRB is FRB 180916 
[81], the host galaxy of which is dramatically different from 
that of FRB 121102. The stellar mass of this galaxy is about 
2 x 10° Mo and the star formation rate (SFR) of the birth re- 
gion is about 0.06 Mayr !. The vast difference between these 
two galaxies suggests that the host galaxies of repeating FRB- 
s span a large range. FRB 190711 has also been localized and 
its host galaxy is also distinct [82]. Recently, it has been i- 
dentified as a repeating FRB [83]. 

The localization of non-repeating FRBs is a big challenge. 
FRB 180924 was localized to a luminous galaxy at redshift 
0.3212 [14]. Using the Deep Synoptic Array ten-antenna pro- 
totype (DSA-10), Ravi et al. (2019) reported the localization 
of FRB 190523 [84]. These two host galaxies are similar to 
each other. The stellar mass is about 10!° — 10!! Mg and the 
SFR is about < 2Mayr !. Several non-repeating FRBs has 
also been localized [15, 82, 85-88]. We list the properties of 
host galaxies in Table 1 and more information can be found 
in FRB host database ". 

The properties of host galaxy can give clues to the pro- 
genitors of FRBs. For example, the host galaxy of FRB 
121102 is similar to that of superluminous supernovae and 
long gamma-ray bursts (GRBs), which has led to the hy- 
pothesis that FRBs are produced by young active magne- 
tars [89-95]. The magnetars are mainly formed through the 
core-collapses of massive stars that can be simultaneously ob- 
served by CHIME and third observing run of the Advanced 
LIGO/Virgo gravitational-wave detectors [96]. Meanwhile, 
FRBs may be produced by pre-merger NS-NS interactions 


It was localized to a low- 
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[97] or activities of newborn magnetars from binary neutron 
star (BNS) mergers [98,99], motivated by the large offsets of 
FRB 180924 and FRB 190523. The offset distribution of B- 
NS mergers derived from new BSE code is consistent with 
the offsets of FRB 180924, FRB 190523 and short GRB- 
s [98]. Besides rapidly spinning magnetars, the connection 
of FRB 200428 with SGR 1935+2154 demonstrates that reg- 
ular magnetars can produce FRBs. Therefore, the FRB host 
galaxy could be Milky Way-type with moderate SFR. Recen- 
t works showed that both intermediate hosts and old hosts 
are consistent with the majority of the FRB hosts [100-102]. 
There is no clear sign of a bimodal distribution in SFR. 


Chatterjee et al. (2017) Spilter et al. (2016) 
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Figure 2 The cumulative energy distributions of FRB 121102 for different 
observations. The value of og is in a narrow range 1.6 < œg < 1.8. (Adapted 
from [103].) 


3 Statistical properties of FRBs 


As the number of FRBs increases, the statistical study of 
observational properties becomes important, especially for 
the classification and comparison with model predictions. 
Many works have been carried out to study FRBs statistical- 
ly, which can give important constraints on the progenitors 
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FRB name redshift | DM (pc em?) | Stellar Mass (x108 Mo) | SFR (Mo yr!) | offset (kpc) | References 
FRB 121102 | 0.1927 | 557 1.4 + 0.7 0.15 + 0.04 0.6 + 0.3 [7] 
Repeating FRB FRB 180916 | 0.0337 | 348.76 21.5 3.3 0.06 + 0.02 5.5+0.1 [81] 
FRB 190711 | 0.5220 | 593.1 8.12 2.9 0.42 + 0.12 3.22.1 [82] 
FRB 180924 | 0.3212 | 362.4 132+51 0.88 + 0.26 3.4+0.5 [14] 
FRB 181112 | 0.4755 | 589.0 39.8 + 20.2 0.37 + 0.11 1.7 x 19.2 [85] 
FRB 190102 | 0.2912 | 364.5 33.9 + 10.2 0.86 + 0.26 2.0 + 2.2 [82] 
FRB 190523 | 0.6600 | 760.8 612 + 401 < 0.09 27 +23 [84] 
Non-Repeating FRBs FRB 190608 | 0.1178 | 338.7 116 + 28 0.69 + 0.01 6.6 + 0.6 [82] 
FRB 190611 | 0.3778 | 321.4 ~8 0.27 + 0.08 11+4 [82] 
FRB 191001 | 0.2340 | 507.9 464 + 188 8.06 + 2.42 1l+1 [13] 
FRB 190614 | ~ 0.6 959.2 ~ 10 [88] 
FRB 190714 | 0.2365 | 504.1 149+ 71 0.65 + 0.20 1.9 x 1.1 [87] 
FRB 200430 | 0.1600 | 380.0 13+6 ~ 0.2 3.0 + 2.4 [87] 


Table 1 Properties of FRB host galaxies. 


[92, 104-111]. It has been found that the statistics of FRBs 
are similar to those of magnetar bursts [106, 111] and solar 
radio bursts [112]. This similarity also supports the associa- 
tion of magnetar bursts and FRBs. 


3.1 Energy distribution 
From observation, burst energy can be calculated through 


B And? F Ve 3 
scc e 
where d; is luminosity distance, F is the burst fluence. Zhang 
(2018) suggested that it should be the central frequency ve 
here instead of the bandwidth of the receiver Av [113], which 
is more reasonable since the FRB spectrum is unknown and 
the emission could extend outside the range of Av. 
Wang & Yu (2017) first found that the energies of 17 burst- 
s of FRB 121102 show a power-law distribution with an in- 
dex of —1.80 + 0.15 [106]. Similar energy distribution of 
dN/dE œ E~! from the Very Large Array (VLA), Arecibo, 
and Green Bank Telescope (GBT) bursts are found by Law 
et al. (2017) [114]. However, a much steeper index of -2.8 
is obtained by Gourdji et al. (2019) using 41 bursts observed 
by Arecibo [115]. Further, a universal energy distribution is 
found by Wang & Zhang (2019) [103]. They collected the 
bursts of FRB 121102 observed by different telescopes, in- 
cluding VLA at 3 GHz [7], Arecibo at 1.4 GHz [5], Areci- 
bo at 4.5 GHz [64], GBT at 4-8 GHz [116], GBT at 2 GHz 
[6, 117] and Arecibo at 1.4 GHz [115]. The threshold power- 
law distribution with a high-energy cutoff is adopted to fit the 
cumulative distribution, which is 
N(» E) = A(El-"* — Ej 2t), 


max 


(27) 


where Emax is the maximumal energy and oz is the index of 
differential distribution dN/dE œ E ?*. The fitting results 
for different samples are shown in Figure 2. It is interest- 
ing that the value of og is in a narrow range of 1.6 — 1.8 for 
bursts at different frequencies, which supports a universal en- 
ergy distribution for FRB 121102. Moreover, similar œg for 
non-repeating FRBs observed by Parkes and ASKAP is also 
found [110,118]. This energy distribution is also similar to 
those of magnetar bursts [106, 111, 119, 120], and solar type 
III radio bursts [121]. A few works have discussed the en- 
ergy distribution across different sources [108, 122-125], and 
the detailed method is introduced in Section 4.1. 


3.2 Width distribution 


Several studies have been performed on the distribution of 
observed pulse widths. Wang & Yu (2017) found that pulse 
widths of 17 bursts of FRB 121102 show a power-law dis- 
tribution with an index of —1.95 + 0.32 [106]. Using more 
bursts of FRB 121102 observed by GBT, a shallower index 
of —1.57 + 0.13 is found [111], which is shown in Figure 3. 
The width distribution of FRB 121102 is consistent with that 
of magnetar bursts [106,111,119]. 
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Figure 3 Cumulative distribution of pulse width for 93 bursts of FRB 


121102. The fit is shown as red line. (Adapted from [111].) 


3.3 Waiting time distribution 


The waiting time At, is defined as the difference between 
occurrence times for two adjacent bursts. It can give the in- 
formation of the burst rate and period. If the burst rate is 
constant, the occurrence rate of waiting time should follow 
the Poisson distribution [126] 


P(Aty) = Ae ^^^», (28) 


where 4 is the burst rate. If the burst rate is time varying, 
the waiting time distribution can be treated as a combination 
of piecewise constant Poisson processes. Generally, for most 
forms of A(t), the waiting time distribution shows a power- 
law form. Wang & Yu (2017) first found that the waiting 
time distribution of FRB 121102 shows a power-law form, 
not a Poissonian distribution [106]. Subsequent work by Op- 
permann et al. (2018) confirmed this result [127]. Figure 
4 gives the cumulative distribution of waiting times of FR- 
B 121102 [111], which is similar to that of magnetar bursts 
[106,111,119] and solar type III radio bursts [121]. It shows 
power-law-like distributions with a slope close to —2 for large 
waiting times. If the variability of burst rate shows spikes like 
ó—functions, the occurrence rate 4 is exponentially growing 
and fulfills the normalization requirement Ld Af (Ada = Ao, 


f(a) =A exp [- ij . (29) 


The waiting time distribution in a given time interval can be 
written as 


Ao 


PM e 
(Aty) (1 + AAt) 


(30) 


The fit of red line in Figure 4 gives Ay = 1.23*050 x 10? ms! 
[111]. 
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Figure 4 Differential distribution of waiting times for 93 bursts of FRB 
121102. The fit is shown as red line. (Adapted from [111].) 


4 Population Study of FRBs 


4.1 Global Statistical Properties 


The all-sky rate is very high for FRBs, typically within the 
range of 10? ~ 10*sky~'day~!. Different surveys give d- 
ifferent rates since their sensitivities and central frequencies 
vary [28]. Basically, the event rate density is more important 
since it is directly related to the source system. It depends 
strongly on FRB luminosity. A recent study suggests that the 
event rate density above 10? erg s^! is ~ 3.5x 10* Gpc? yr^! 
[122]. The precise cosmological volumetric rate is waiting to 
be revealed since the FRB redshift distribution V(z) is still 
largely unknown. Except for 13 FRBs that have been local- 
ized with identified host galaxies [9, 14, 81, 82, 84-87], the 
majority of FRBs are lacking in direct distance information. 
Several pioneering works have put constraints on V(z) with 
current statistical properties of FRBs [107, 128-131]. 

The event rate theory is briefly introduced as follows. Oth- 
er than V(z), the frequency-dependent FRB detection rate 
is largely determined by the intrinsic luminosity function. 
We can define a function O(L,,z) to incorporate both fac- 
tors [107]. The form of O(L,,z) is unclear from the first 
principle, however we can model it with some simplified 
assumptions. If we assume that the luminosity function 
does not evolve with redshift, O(L,, z) is therefore separa- 
ble O(L,, z) = O(L,)¥(z). Further, different models could be 
applied to ®(L,) and V(z) respectively. For instance, sim- 
ple luminosity function models like standard candle, power- 
law function and log-normal distribution have been adopted 
in previous studies [129, 131]. More specially, the Schechter 
function has been discussed [108, 132]. For the redshift evo- 
lution, the most commonly adopted model is that Y((z) traces 
the cosmic SFR [129, 131]. Note that there might be a time 
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delay of ‘¥(z) from SFR if an FRB is produced by compact 
binary mergers [107, 118, 133]. This delay effect has been 
taken into account for short GRB before [134, 135]. Also, 
other models such as a constant comoving density or ¥(z) in 
proportion to cosmic stellar mass density have been discussed 
[131]. 

The generic mathematical formalism that relates the ex- 
pression of O(L,,z) with observational properties has been 
completed in Macquart & Ekers (2018b) [107]. Usually, the 
cumulative distribution of FRB flux density/fluences (known 
as log N — log S distribution) is devoted to constrain G(L,, z). 
If most of FRBs are in our local universe, the slope of 
log N — logS plot should be —1.5 [136]. However, differ- 
ent surveys give different power-law indexes in a wide range 
[129, 130, 137-140]. Meanwhile, the cosmological (V/Vimax) 
tests of different samples deviate from 0.5, which disfavor 
an Euclidean geometric distribution [133, 141]. Furthermore, 
the DM distribution of FRBs can constrain G(L,, z) due to the 
relation between DM and z in Eq. (7). Assuming the models 
of V(z), several works have constrained the FRB luminosi- 
ty function, or equivalently, the enengy distribution function 
[108, 110, 118, 123, 128, 131, 142, 143]. Recently a Bayesian 
method has been applied to give both the best-fitting val- 
ues and error-bars of the parameters in Schechter luminosity 
function [110, 122]. Besides, the intrinsic spectra of FRBs 
can influence the detection rate, since the emission frequency 
in source frame differs from that in the observer frame due 
to cosmological expansion. This effect is embedded in a k- 
correction factor (1--z) * if the power-law spectrum L, œ y ? 
is assumed [107,113, 131]. There are a few works that have 
discussed the impact of index «œ [128, 131, 133, 143]. Howev- 
er, currently it is poorly constrained from observations except 
for some small samples [144]. 


4.2 Repeating and Non-repeating: A Single Population? 


There are hundreds of FRBs from different surveys now 
recorded by the Transient Name Server ?, among which most 
are one-off events and only twenty are identified as repeater- 
s. Generally, an FRB is recognized as a repeater if a second 
burst has been detected in the same position of the sky, which 
relies much on the precision of localization. However, there 
are several factors that make a repeater look non-repeating. 
First one is the selection bias of the receiver. If the pulse 
width of the second burst is close to the instrumental sam- 
pling or smearing timescale, then this signal is probably be- 
ing missed [145,146]. Also, the beam shape could play a role 
[147]. Second, the energy distribution function of individual 
FRBs matters a lot. If the second burst is much dimmer than 
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the original one, it may fall below the detection threshold. 
This has been implied by the repetition of FRB 171019, for 
which two fainter bursts by a factor of ~ 590 than the origi- 
nal burst that have been detected [148]. It may be just a lucky 
case and we could probably have missed many faint bursts of 
other FRBs. Third, the repetition statistics could differ a lot 
for each repeater. The waiting times for the majority of FRBs 
may be much longer than those of the frequent repeaters like 
FRB 121102. This problem will be alleviated as the observ- 
ing time accumulates. 


These above lead to a debate: do all FRBs repeat [145, 
149, 150]? There is no firm conclusion but some emerging 
clues seem disfavoring a single population. The burst mor- 
phology looks different for repeaters and non-repeaters. Re- 
peaters like FRB 121102 have wider pulse width and show 
complex sub-burst structures [11-13,56,151], while pulses of 
non-repeaters are usually narrower [13, 143]. If the temporal 
and frequency structure is caused by self-modulation during 
propagation [152], the circumburst environment of repeaters 
and non-repeaters should be very different. Further, the po- 
larization properties seem different, though not prominent in 
consideration of the presence of FRB 180301 [69]. Moreover, 
there are some differences in the mass and SFR of their host 
galaxies, however, a clear sign of dichotomy has not been 
found [15, 100]. Note that these clues suffer from a small 
number of repeaters. With more repeaters coming soon, bi- 
modal distributions are expected in these statistical properties 
if repeaters and non-repeaters are physically different. 


Still, it is fair to argue that all FRBs repeat on the suppo- 
sition that their repeating modes are diverse. Several work- 
s have tried to reproduce the observed statistical properties 
with different models of repeating [153, 154]. For instance, 
the difference in pulse width distribution can be ascribed to 
FRB beaming [155]. Two factors are decisive in repetition 
modeling: the burst energy function and the distribution of 
waiting time [123, 147, 150, 156]. Further, assuming the cos- 
mic evolution of comoving FRB number density, we can cal- 
culate how many repeaters will be identified for a given sur- 
vey. The energy function is commonly assumed as power- 
law or Schechter forms [123, 147, 150, 156]. Since the dis- 
tribution of waiting time of FRB 121102 is non-Poissonian 
[18, 111, 115, 127, 157], a Weibull distribution is usually 
adopted [123, 127, 147, 150,156]. A lack of repeaters in sev- 
eral surveys suggest that most FRBs should not repeat as fre- 
quently as FRB 121102 [145, 158, 159]. Generally, with the 
accumulation of the observing time, the observed repeating 
fraction (defined as the number ratio of repeaters among all 
FRBs) is expected to increase to nearly 100% [123,150,156]. 
However, if there exist physically non-repeating FRBs, this 
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fraction will peak at a certain observing time and the peak 
fraction is less than 100% [153, 156]. The current total ob- 
serving time of CHIME may be still insufficient to bring out 
the peak. Also, the repeating fraction will decrease with the 
increasing redshift since the repetitions of FRBs with closer 
distances are easier to be detected [123,153]. This means that 
the DM distribution of repeaters should concentrate on a low- 
er value than non-repeaters. Observationally, this trend has 
not been confirmed yet since the number of repeaters is lim- 
ited. A bottom line is that the volumetric rate of apparently 
non-repeating FRBs exceeds those of candidate cataclysmic 
progenitor events [160], therefore, there are many underlying 
repeaters waiting to be discovered in the near future. How- 
ever, the existence of a physically one-off population can be 
confirmed once the direct association of an FRB with a cata- 
clysmic event is observed. 


4.3 Periodicity 


The periodicity of FRBs is a new character discovered in ear- 
ly 2020. The repeat bursts of FRB 180916 arrive regularly 
with a period of 16.35+0.15 days [16]. More intriguingly, the 
activity window is found to be narrower and earlier at higher 
frequencies recently [78]. Further, a period of ~ 157 days for 
FRB 121102 has been claimed by two different groups using 
different burst data [17,18]. Currently, period finding has on- 
ly been performed to three relatively “bursty” repeaters (the 
above two together with FRB 180814), however, no signif- 
icant periodicity has been found for FRB 180814 yet [161]. 
Usually at least tens of repeat bursts are needed to achieve 
high enough confidence level after period folding process. 
This implies that there may be underlying periodicity for oth- 
er repeaters. The situation is very much analogous to that of 
repetition two years before. Future debate might arise like 
this: do all repeaters have a periodicity? One thing is obvi- 
ous that the unique period is distinct for each individual case, 
which may vary from seconds to years. Learning from repe- 
tition, we can do period modeling by assuming a distribution 
of period, which is left for future work. 

Theorists have proposed several possible physical scenar- 
ios for the periodical activity [162]. Generally the existing 
models can be classified into two categories. First kind of 
models attributes the periodicity to the orbital motion of a 
binary system [97, 163-170]. This binary usually involves a 
neutron star (NS) as FRB emitter together with a compan- 
ion. The companion can be an O/B type star [166, 168] or 
a compact object [97, 164, 165]. There are multiple ways to 
achieve periodicity and one possibility was proposed that the 
wind of the companion can pave a way for the FRBs to es- 
cape [165,166]. This opened window would point to us pe- 
riodically due to the orbital motion. However, this scenario 
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has recently been disfavored since it predicts a wider activity 
window for higher frequencies [78]. Alternatively, there may 
be an asteroid belt surrounding the companion and the NS 
collide with them in a certain orbital phase [163, 171]. The 
second kind of models suggests that the periodicity is due to 
NS precession [172-177]. The NS could either undergo a free 
precession [172, 174, 175] or forced precession [173], while 
the latter may be caused by the surrounding disk [176, 177]. 
However, one important issue for precession models should 
be addressed properly: why are FRBs produced at a fixed re- 
gion on the NS surface? Instead of an NS, the precession of 
a jet produced by a massive black hole has been discussed 
[178]. Note that a third option that associates the periodici- 
ty with the rotation of a magnetar has been proposed [179], 
however, there is no evidence that such a slowly-rotating 
magnetar exists and is able to produce FRBs. The period- 
ical behavior is closely related to the radiation mechanism 
and source model of FRBs, which will be discussed in next 
section. 


5 Radiation Mechanism for FRBs 


The radiation mechanism is always the core issue of a new 
astrophysical signal. It may take a while before we can figure 
out the mechanism of FRBs. However, we can learn a les- 
son from similar phenomena like pulsar radio emission and 
GRB. Very recently, Zhang (2020) gives a very nice review 
on this issue [30]. The high brightness temperature of FRBs 
from Eq. (2) implies that the radiation must be coherent. Ba- 
sically there are three main approaches to generate coherence 
(for a review, see [180]). The first one is coherent curvature 
emission by bunched particles (also known as antenna mech- 
anism) [95, 181-183]. The second mechanism is relativistic 
plasma emission by reactive instabilities (also called plasma 
maser in some literatures). The particles' kinetic energy is 
transferred to plasma waves through a beam instability and 
is finally converted to escaping radio emission [184, 185]. 
The third mechanism is maser (microwave amplification by 
stimulated emission of radiation). The radiation can be am- 
plified if “negative absorption" is reached. The key point is 
that the number of high-energy electrons should exceed that 
of low-energy ones under certain physical conditions, which 
is called “population inversion”. Till now various versions 
of maser has been proposed for FRBs, including plasma syn- 
chrotron maser, vacuum maser and synchrotron maser from 
magnetized shocks [59, 90, 186-188]. 


5.1 Coherent Curvature Emission 
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5.1.1 Synchro-curvature emission from a single charge 


In vacuum, the trajectory of a relativistic electron in the pres- 
ence of a magnetic field can be decomposed into two com- 
ponents: a motion along the field line and a gyration in the 
perpendicular plane. The latter leads to synchrotron emis- 
sion and the former could produce curvature emission if the 
field line is curved (e.g., a dipole field geometry). Therefore, 
the realistic emission of relativistic electrons in a curved field 
can be generalized as synchro-curvature radiation, which is 
illustrated in Fig 5. 


Figure 5 Schematic picture of synchro-curvature radiation. 


Suppose the magnetic field is along $-direction B = Boe; 
and its curvature radius is p. A relativistic electron moves 
with a Lorentz factor y and the angle between v and B is a. 
The guiding center drifts along the field line with an angular 
velocity Qo. The gyro-radius is defined as 


yv,m,c yp sin amec? 


eBo = eBo 


-csina/op, (31) 


where wg = eBo/(ym,c) is Larmor frequency. 
The energy radiated per unit solid angle per unit frequency 
interval of a moving charge is [189] 


2 


f nx (nx peda. (32) 


dE eo 


dwdQ — 472c 


We can construct a Cartesian coordinate system, making sure 
that x—y plane overlaps with p—¢ plane and initially the guid- 
ing center is on the y-axis (Fig 5). Without loss of generality, 
the unit vector pointing to the observer n is in x — z plane 
and makes an angle 0 with the x-axis. Substituting the ex- 
pressions of n, B, r in the xyz coordinate and integrating the 
solid angle, we can obtain the characteristic frequency and 
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radiation spectrum of synchro-curvature radiation [190, 191], 
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where rž ~ c?°/[rgwh + (p + rg)Q2], and 


2 2 
ry + pre —3 3 
05 = (45 cos? a eos + — cos œ cos 0 
P P 
MEE 1 
+— sin asin) —. (34) 
TB TB 


The above formulae can be easily degraded to synchrotron 
and curvature radiation case. For pure synchrotron radiation, 
p = co, a + 0, Q = 0, Q» = sin’ a/rp, and r2 Q2 = 1. So, 


3 : 
We 2* Wp SIN Q, 
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Then for pure curvature radiation, a = 0, rg = 0, Qo = c/p, 
Q» = 1/p, and r? 03 = 1. So we have 


3 
We = = 
a” 
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If FRBs are produced inside the magnetosphere of NS 
where B > 10°G, the gyration damps very quick since the 
synchrotron cooling timescale is proportional to B-*. The 
electron moves along the field line and only curvature emis- 
sion is important [192]. However, outside in the low B region, 
both curvature and synchrotron emission can be effective and 
Eq. (33) should be applied. 


5.1.2. Coherent curvature emission by bunches 


Now we can discuss the radiation spectrum from a bunch of 
N particles in the magnetosphere. Coherence can be achieved 
if the phases of EM waves emitted by each individual electron 
are near the same. According to Eq. (32), the total radiation 
intensity of N electrons is 


N 2 


d Eit E eo? C p 
dwdQ  4mc |J ; 


n x (n x Bye" "Ot ,(37) 


where the subscript j represents the j-th electron. This inte- 
gration can be simplified under certain assumptions and has 
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been discussed in detail by Yang & Zhang (2018) [95]. Here 
we list some of their important results below. 

Considering a one-dimensional (1-D) bunch with finite 
length L moving along a magnetic field line of curvature ra- 
dius p, the position vector from the j-th electron to the first 
one is Ar,(t) = rj(f) — r(t). Since the emission is highly 
beamed, this vector could be considered as time-independent. 
Eq. (37) can be rewritten as 


2 
d Eit eo? 


+00 
[ nx (n x Byet-nrolog, 


dwdQ ^ 472c 
" 2 
x ee) (38) 
J 
The phase stacking term F,, = ba ew? can be e- 


valuated analytically. Defining a typical frequency wz 
2c/(L cos 0) where 0 is the observing angle, we have [95] 


N?, W K wz, 
ns MM (39) 
N (2) , WL KWX Weoh- 


where Weoh ~ (ILP wL is a frequency boundary beyond 
which the emission becomes incoherent. Further, if we as- 
sume the electrons have a power-law energy distribution, i.e., 
Ny) x y? where ym < y < Ymax, the resulting radiation 
spectrum has a break at wm = We(Ym) given by Eq. (36). 
The power-law index of the spectrum turns from 2/3 (for 
w < Wm) to —(2p — 4)/3 (for w > wm). For more compli- 
cated scenario, a 3-D bunch with a half opening angle y has 
been considered. If y is larger than the angle of beamed cone, 
only part of the radiation by this bunch can be observed. This 
also leads to a spectral break from 2/3 to 0, at the frequency 
Wy = 3cK(pg?). 

Taking all these above into account, the total radiation 
spectrum of a 3-D bunch filled with power-law distributed 
electron has a four-segments-broken-power-law shape, char- 
acterized by three frequencies wz, Wm, Wy. The indices of 
four segments change according to the relative values of 
WL, Wm, Wy and the complete analytical expressions can be 
found in Yang & Zhang (2018) [95]. 


5.1.3 Bunch formation 


Coherent curvature emission by bunches has long been pro- 
posed as a possible mechanism for pulsar radio emission 
[193-195]. However, it remains controversial how the bunch- 
es could form effectively, and consequently some critiques 
has been put on this mechanism [196-198]. It seems unlike- 
ly to solve this leftover problem completely in FRB field but 
some possible ways could be explored. 
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For instance, the possibility of bunch formation by two 
stream instability under the physical condition of FRBs has 
been discussed [181, 188]. Instability caused by counter- 
streaming electrons and positrons leads to density fluctua- 
tions. Assuming that electrons and positrons have the same 
density and the four velocity distribution f(u) is symmetric, 
the dispersion relation can be written as [188] 


+00 442 +1 
fw = 2 FP) = 
borm hor em 


where k is wave-number and 8 = u/y is velocity. Further, as- 
suming a step-function distribution f(u) within velocity range 
Bmin < B < Bmax, the instability growth rate Im(w) can be 
determined, which is given by the imaginary part of com- 
plex frequency w. Defining an effective plasma frequency 
Wpett = Gp (y 2)? where (y?) = [7^ y f(u)du, the effec- 
tive skin depth is then fskin = C/Wp er. Large density fluctua- 
tions can be produced on length-scales larger than fskin, 1.e., 
bunches with L > fskin could possibly form via two stream 
instability. Note that above treatment is greatly simplified s- 
ince f(u) is highly uncertain in reality. More future work is 
needed to address the issue of bunch formation. 


5.2 Relativistic Plasma Emission (Plasma Maser) 


The theory of plasma emission was first proposed to explain 
solar radio bursts [199]. Later its relativistic version was ap- 
plied to pulsar radio emission [200]. Generally plasma emis- 
sion process includes two stages: Langmuir-like waves are 
generated by a beam instability and then part of the wave en- 
ergy is converted to escaping radiation [184]. 

The beam instability could occur when a beam of relativis- 
tic particles runs into a background plasma. The growth rate 
of beam instability depends on the dispersive properties of 
the plasma. To ensure effective plasma emission, the growth 
rate should be relatively fast. Assuming the four-velocity 
distributions of the beam f,(up) and target plasma f(u), the 
dielectric tensor of the plasma can be obtained and the dis- 
persion relation of Langmuir-like waves can be determined 
[184]. The growth rate is enhanced at the resonant frequen- 
cies. The Cherenkov resonance condition reads 


1 - nß cos 0 = O, (41) 


where n is the refraction index. Also, the cyclotron- 


Cherenkov (or anomalous Doppler) resonance condition is 
w(1 — nf cos 8) + cg = 0. (42) 


Plasma maser is effective if the growth rate at resonant con- 
ditions is fast enough. However, there are two main diffi- 
culties for this mechanism. First, it remains questionable 
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whether Langmuir-like waves through beam instabilities can 
be generated in the NS magnetosphere [185]. Second, usu- 
ally the growth rates for beam instabilities seem too small 
to allow effective wave growth [180]. In the case of FRB 
121102, cyclotron-Cherenkov resonance condition requires 
unrealistically high beam Lorentz factors and the growth rate 
of Cherenkov instaiblity is too small [188]. Thus, plasma 
maser can hardly explain the emission of FRB 121102. For 
completeness, however, we list this mechanism here since it 
represents the second way of generating coherence. The liv- 
ing space is limited even it might be applicable to some par- 
ticular low-energy bursts. 


5.3 Plasma Synchrotron Maser 


Synchrotron self-absorption can happen when the emitted 
synchrotron radio waves propagate in a weakly magnetized 
relativistic plasma. The self-absorption coefficient can be ex- 
pressed as [201, 202] 


c EE 


E=- gm 0 dE | E? | dv 
&c (* N(E) d|E?dP/do)] 
- —— dE 4 
w Jo E dE , (43) 


where E = ym,c? is electron energy. Normally u is positive 
but it can turns negative if the plasma is dense and relativistic. 
A direct requirement is that + E | > 0, i.e., the energy dis- 
tribution of electrons grows faster than E?. This means that 
population inversion is necessary. Further we need to explore 
at which frequencies this negative absorption is significant. 
The synchrotron radiation power in the presence of a back- 


ground plasma is [203]. 


dP 3e | i 
_ V3e yop sina w f Baoi (44) 
[w 


d 3m 
2 2ncv 1 Ry s 


where the critical frequency o. is then 


, 


3 
w, = 5798 sina(1 + y we [ap ?P., (45) 
If y’w;/w* < 1, Eq. (44)(45) tum back to Eq. (35). On the 
contrary, the influence of relativistic dense plasma (y, wp are 


large) should be taken into account. Substituting Eq. (44)(45) 
into Eq. (43), we get [204] 


d 
po ET [E?dP/dw] œ z ? dX), (46) 
where 
_@ 2 w 
iia w.  3o?opsina 
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Ëz) = 2z j Ks/3(x)dx — 2°Ks/3(2). (47) 


The function ®(z) has a negative minimum value of ~ —0.24 
at the frequency vy: = vp(vp/ vg)! ? [205]. For v < vg, u is 
always negative. 

The intensity /,, after radiation transferring in a source of 
size L is 


I, = Joc -e"), (48) 
H 


where j, = T (dP/dw)N(E)dE is the total emissivity and 
T = pL is the optical depth. As long as u is negative and 
lulL >> 1, the radiation intensity is amplified exponentially 
[204] 


I = Gol|ul) exp(ulL). (49) 


This indicates that a significant fraction of electron energy 
is converted to the strong synchrotron maser emission. Note 
that |u|L > 1 is achieved only for a narrow range of frequen- 
cies, so the maser emission has a narrow spectra [204, 205]. 
Also, its typical frequency is lower than synchrotron self- 
absorption frequency [204]. In some circumstances, this plas- 
ma synchrotron maser could arise around GHz, thus being 
regarded as a possible mechanism for FRBs [186]. 


5.4 Vacuum Maser 


For synchrotron emission in vacuum, if Eq. (35) is substi- 
tuted into Eq. (43), we always have arabian > 0 and thus 
u > 0. However, we should notice that Eq. (35) has been in- 
tegrated over solid angle. Actually, the synchrotron radiation 
power is angular-dependent and iV is relevant. At some 
certain observing angle 0, the self-absorption coefficient can 
be negative [187,188, 206]. 


The emissivity for a single electron is [207] 
dP 9ercB? ( v V 2,2 
— = —] a 0 
dvdO ^ 64m vg (=) harass 
[+ 7°0)K330) + € K?40)]. 60) 
where vg = wg/27, v, = wc/27 and y = 5-1 +y@)*, The 
angular-dependent absorption coefficient can be rewritten as 
1 ? N(y) à : 
Z L ypj)dy, — 61) 
yp oy 


where p = (y? — 1)!/? is the electron momentum in units of 


Jy) = 


u(0) = 


2mev? Jy 


mec. A differential absorption cross section has been defined 
as [206] 


do(0) L lð 
- 2 
dd GaP op Oy Lypjv()] (52) 


Substituting the expression of j,(0) and letting t = y6, 
do,(@) 2e 1 
dQ — 9 By sin? 


— (13 - LD + 0K30) 
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+t(11t — 1)KÎ 30) - 6y(¢ - 2) 
x[C1 + 1)K2/3(y)Ks/3(y)  tK1750)K4730)1)- (53) 


This cross section is negative only if both à > V2/y and 
v/ve > (yo)? are satisfied [187,206]. The maser emission 
can be observed outside the 1/y emission cone of the rela- 
tivistic electrons. For lower frequencies v/v, « (y0)? or in 
the 1/y cone the cross section turns positive with a larger ab- 
solute value [187]. This implies that the distribution of pitch 
angles should be narrower than 1/y, otherwise any photon 
will be preferentially absorbed by a particular electron whose 
emission cone wraps the photon's moving direction. Howev- 
er, it is unclear whether such narrow pitch angle distribution 
can be achieved. 

Other than vacuum synchrotron maser, the vacuum curva- 
ture maser inside the magnetosphere has been considered. S- 
ince the expressions of emissivity for synchrotron and curva- 
ture radiation are very similar, the above procedure of calcu- 
lation has been repeated [188]. The curvature self-absorption 
cross section could be negative for 0 < Qy))!. However, 
since electrons move along curved field lines in the magne- 
tosphere, the curvature photons emitted by an electron will 
be inevitably absorbed by electrons on other field lines since 
their intersection angles are large. Therefore, it is unlikely for 
vacuum curvature maser to be responsible for FRBs. 


5.5 Synchrotron Maser from Magnetized Shocks 


Relativistic perpendicular magnetized shocks has been stud- 
ied extensively in literatures. Traditionally, a set of conserva- 
tion equations of physical quantities between upstream and 
downstream plasma can be written. However, these macro- 
scopic conservation laws fall short of describing the mi- 
croscopic behavior of the plasma, such as the distribution- 
s of particle momentum, energy partition between different 
species and so on. These kinetic issues shall be addressed by 
collisionless plasma physics. The development of particle- 
in-cell (PIC) simulations since early 1990s greatly helps un- 
derstand the plasma properties. In principle, particles reflect- 
ed at magnetized shock front will achieve bunching in gyra- 
tion phase and synchrotron maser (or relativistic cyclotron) 
instability will develop [208]. Strong coherent EM signal can 
be produced and the reabsorption of these cyclotron waves 
serves as the main heating mechanism of downstream plas- 
ma [209]. This EM signal occurs before the radiation from 
heated downstream particles so that it is called a precursor. 
The FRB emission is probably just this kind of EM precursor 
[90]. Here we will revisit the results of 1-D PIC simulations 
and briefly introduce the underlying physics. 

First we present the analytical expressions of shock jump 
condition for an ideal MHD plasma. The conservation of 
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particle number, energy, momentum and magnetic flux are 
respectively expressed as [210] 


N2 1+ Pshock 


z= mmh, 54) 

Ni Pshock 
B 

(1 + Bsnoa) Numyic (1. 0) = Bsnock [e 5 =] j^ 55) 
B 

(1 + Baa) Nimyyc (1 +0) = pa + = (56) 
B 1 snoc. 

Bi _ 1+ Bsnock (57) 


B 1 P shock 


where subscripts 1, 2 represents upstream and downstream 
respectively. The upstream is assumed to be cold and rela- 
tivistic, p; = 0, B1 = 1. For the hot relativistic downstream 
plasma, e? = p» KE — 1) where Î is adiabatic index. The 
nav : B 
magnetization parameter is defined as o = ;——15—. From 
zaNumc^yi 
above equations we can obtain the equation for shock veloc- 


ity shock: 
c 


1 
Pd |. a -l= 
( + - Jas 2 


which has a solution, 


1 {(5 ==) 
+ + 
20 + 1/o) | \2 c 


r 11^ 


The ratio of the shock jump is 


^ 


pie 1] Ass : ( : 5) = 0,(58) 
o 2 


PB shock = 


N B 1 
Sui c4 x . (60) 
Ni Bi P'shock 
kT. 1 + Bshoc 
ae reor e] (61) 
myic 2B hock 


where T» is downstream temperature. Substituting Eq. (59) 
we can express this ratio as a function of c and f. 

However, the results from PIC simulations shows that both 
No/N, and KT; /my\c? deviate from Eq. (60)(61). The reason 
is that the ideal MHD formalism does not account for the ad- 
ditional wave fluctuations that can dissipate the flow energy. 
The upstream field fluctuation is due to the EM precursor, 
which can be expressed with a parameter £, 


WOW E. Ex 
—(6E;, -óB?) x ——(óEj,0Bj) = —£&B, 62 
g, (El + Bi) ~ —7-(6E10B1) = By, — (62) 


Besides, there are two sub-dominant fluctuations in the down- 
stream. The parameters ¢, 77 are introduced to describe its EM 
and electrostatic fluctuations respectively, 


1 1 1 1 
zz OE + 6B3) = a g OE) = iB: (63) 


Taking these modifications into account, Eq. (55) and (56) 
can be rewritten as 
1 
Nymyie*(1 + Bao) + z= BT + snoek) = EC — soa] 
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1 1 
- je reaps e+ 7)Bi| Bshock, (64) 


1 1 
Nimy ic? + — Bi qd + Bshock) + —éEBi(1 — Bshock) 
4n 4r 
1 1 
= — B2 + —(ć -9m. 
Pat 5B, + gé Br») 


while Eq. (54) and (57) remain unchanged. Then the equa- 
tion for Bshock becomes 


1 " a fl 1 ^ | 92 
( T F - £) Bee + a-b[; T c e£ dety Pstoa 


^ 


2-T 
9 rose B EM = 0.(66) 


1 
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The form of density jump Eq. (60) remains while the down- 
stream temperature is now 


KT» 1 ro = 1 all 


Nn: = Bshock 
my|c 


2 
shock 


op. shock 
1+ Bshock 


Three efficiencies are defined as the fractions of flow energy 
that carried away by these fluctuations, 


l= shoc 
p= 3 Bshock _  & (1- 


[£(1 — Bshock) — £ + n). (67) 


2Nı 
N2 


1+1/71+Bshock 14+1/o 


on PB shock = 6,7 M 
1+1/0 1 + shock 1+1/o0 M 


) . (68) 


Il 


Sin (69) 
Simulation results show that fe, < fe [210], and we are only 
interested in the EM precursor. fg has a peak value around 
10% at o ~ 0.1, and increasing o leads to decreasing fg. For 
very high magnetization o > 1, fg has an asymptotic form 
of 7 x 107^/a? [211]. The spectrum of EM precursor peaks 
at Wpeak ~ 3wp max[1, vc] (in the post-shock frame) and is 
relatively narrow-banded with Aw/Wpeak ~ a few [211]. 

The synchrotron maser instability can be understood as 
follows. Bunching in gyration phase of incoming particles 
has been verified by PIC simulations [208, 210-212]. For 
pure electron-positron plasma, this means e* gyrates in the 
magnetic field with the same Lorentz factor. Therefore, the 
distribution of pairs in momentum space can be described as 
a cold ring, which reads 


1 
fuu) = TER ó(u, — uo)ó(uj), (70) 
Tuo 


where the symbols || and are with respect to the magnet- 
ic field direction. The population inversion has been reached 
since low-energy pairs are absent. If the magnetic field is as- 
sumed on z-axis and cyclotron waves from pairs propagate on 
x-axis, the dispersion relation can be written as 


mA, (71) 
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where e is the dielectric tensor. Sustituting Eq. (70), the four 
components of € can be expressed in analytical forms [209]. 
The dispersion relation shows two unstable branches corre- 
sponding to EM waves and magnetosonic waves [209,212]. 
It has been demonstrated that the growth of EM waves can 
be very effective at low harmonics thus a synchrotron maser 
(EM precursor) can be produced [209, 212]. Simultaneous- 
ly, the energy spectrum of downstream pairs is heated to a 
relativistic Maxwellian distribution [209,210]. Note that the 
above discussion can be generalized to an ion-pair plasma 
though the expression for s becomes more complicated. In 
the presence of ions, a high-energy non-thermal tail of ener- 
gy spectrum of pairs appears, which is due to the absorption 
of the ion cyclotron waves [212, 213]. 


To conclude this section, we should be aware of the fact 
that above mechanisms do not include all, since more than 
a dozen of pulsar radiation models have been proposed and 
developed. It is unclear which mechanism would work for 
FRBs. According to the site where these mechanisms take 
place, they are classified into two types: close-in (i.e., insid- 
e the NS magnetosphere) and far-away (i.e., outside the NS 
magnetosphere) (see the latest review of Zhang (2020) [30]). 
A few works have been devoted to discuss the applicability 
of different mechanisms to individual bursts. As an example, 
a critical analysis on radiation mechanism of FRB 121102 
has been nicely performed by Lu & Kumar (2018) [188]. 
They concluded that given the properties of FRB 121102, 
plasma maser and vacuum curvature maser have been exclud- 
ed. Vacuum synchrotron maser may survive, but needs fine- 
tuning. Plasma synchrotron maser and synchrotron maser 
from magnetized shocks suffer from the problem of low ra- 
diation efficiency. Coherent curvature emission satisfies the 
observations best, however, the relevant magnetic reconnec- 
tion physics is highly uncertain. For other FRBs with lower 
energies, the constraints on radiation mechanism may be to- 
tally different. For recent FRB 200428, models based on both 
coherent curvature emission and synchrotron maser emission 
have been applied successfully [124, 214-216]. However, di- 
verse PA swings of FRB 180301 clearly favor the close-in 
mechanism [69]. One thing that we know from the GRB field 
is that the radiation mechanism is the kernel of FRB physic- 
s and it is not easy to be made clear completely. Similar to 
GRBs, the controversy on radiation mechanism may last a 
while before reaching a final conclusion. A good mechanism 
needs to explain FRB energetics, duration, spectrum, pulse 
structure, polarization, repetition and periodicity, hence plen- 
ty of work should be done in the future. 
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6 Magnetar Source Models 


The origin of FRBs has been an outstanding issue ever since 
its discovery. Various source models have been proposed, 
most of which involve compact objects (see [217] for a the- 
ory catalogue). Different models predict different EM coun- 
terparts and sometimes multi-messenger signals [218]. Note 
some progenitor models are catastrophic and can only be ap- 
plied to one-off FRBs. The recent discovery of FRB 200428 
has confirmed that magnetars are capable of producing FRB- 
s. Moreover, localized FRBs seem to be consistent with mag- 
netar progenitors in host galaxy properties [101, 154]. Al- 
though it is still early to say that all FRBs are from mag- 
netars [102], at this time magnetar models attract the most 
attention of the community. In this section we will introduce 
four well-developed magnetar models and their predictions 
[59, 171, 219, 220]. 


6.1 Afar-away model of Metzger et al. (2017, 2019) 


The model proposed by Metzger et al. is based on syn- 
chrotron maser emission from magnetized shocks, which is 
far away from the central magnetar. The picture is shown in 
Figure 6. A magnetar is embedded in a supernova (SN) envi- 
ronment and the wind blows a bubble (nebula) behind the SN 
ejecta. Matter ejected during magnetar flare is composed of a 
highly relativistic shell and a subrelativistic ion tail. The rela- 
tivistic ejecta shell collides with the leftover ion tail from the 
previous flare and forward-reverse shock system is formed. 
The ion tail is magnetized with assumed ø ~ 0.1 — 1, which 
is beneficial for producing synchrotron maser by the forward 
shock. Assuming the ion tail has a radial density profile 
Next œ r*, the dynamical evolution is then quite similar to 
that of GRB afterglow [221]. 


jabu SN ejecta 


relativistic 
flare 


ion shell 


engine 


z 


=A ~ cót— 


UT twAT—_. Ry Rej 


Figure 6 
[59]). 


The schematic picture of Metzger et al. model. (Adapted from 
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The density of the unshocked ejecta shell is 


[2 1 


"M 72 
ôt Anr?myc?T 2. (a) 


Mej 
where Tej is the initial Lorentz factor of the shell and €, ôt is 
the energy, duration of the flare. The reverse shock crosses 
the shell at a time fa& ~ ôt [221]. For the typical flare pa- 
rameters, the reverse shock is always relativistic. Before the 
crossing, the Lorentz factor of the shocked region is 


1/4 


fl T 
T= 1 ort, (73) 


where f = ng/ney oc r-2 After crossing, the Lorentz fac- 
tor approaches the Blandford-McKee self-similar form T° œ 
7-3? [222]. Substituting t ~ r/(2cT?), the time evolution of 
Dis 


(k-2) 
LOV, t< ot 
pao, t> ot 


As discussed in Section 5.5, the synchrotron maser peaks at 
Wpeak ~ 3Wp max[l, Yo] in the comoving frame. For the 
assumed o < 1, the peak frequency in lab frame is 


1 
Vpk & 54; OTO) (75) 


However, sometimes the observed FRB frequency may be not 
this intrinsic peak frequency. The opacity of induced Comp- 
ton Scattering (ICS) can be very large for the low energy part 
of the maser spectrum. The observed peak frequency is shift- 
ed to a higher value set by Vmax = V(t. = 3). Substituting Eq. 
(74), time evolution of the observed FRB frequency is then 


" 
fH, t < ot 
5/4,1/4 ^ ES > 

vma T va ur MM (76) 
E CWD, tz ot, 


This naturally predicts the down-drifting of the FRB emis- 
sion. Since there is strong ICS absorption, the effective radi- 
ation efficiency is even lower than maser efficiency of ~ 107°, 
typically ~ 1076—1075. Another prediction is that the acceler- 
ated electrons in the downstream can produce a high-energy 
(keV-MeV) counterpart. Its duration is of order ~ millisec- 
onds, and the luminosity and spectrum highly depends on the 
flare parameters. 

This model has been successfully applied to FRB 200428 
in different details [215,216, 223,224]. The observed low ef- 
ficiency Eradio/ Ex ~ 107? of FRB 200428 matches the mod- 
el prediction well. The duration and luminosity of the as- 
sociated XRBs can be explained [215]. Especially, by in- 
troducing a density-jump structure of upstream medium, the 
double-peaked character, luminosity ratio and emission fre- 
quency of two radio pulses can be well explained [216]. 
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Moreover, three X-ray peaks are expected [216], which is 
also consistent with observation [22, 23]. However, obser- 
vationally the X-ray peak energy seems to increase with time 
[23], which needs further consideration. Non-thermal accel- 
eration of electrons can be realized if sufficient ions exist in 
the downstream [212]. The superposition of two components 
could change the peak energy. Besides, the XRBs may sim- 
ply have a different origin from FRB [124, 223, 224]. This 
model also requires a large amount of baryonic matter (about 
0.005 solar mass), which may be larger than the typical mass 
of a magnetar outer crust [224]. Note that this model needs 
upgrade in order to explain the complex PA behaviour of FR- 
B 180301, since the magnetic field configuration is basically 
fixed and a flat PA curve is expected. 


6.2 A far-away model of Beloborodov (2017, 2020) 


The radiation mechanism in Beloborodov model is also syn- 
chrotron maser emission, however a main difference from 
Metzger et al. (2019) is the upstream medium. In this model, 
the emission region is also far away from a central magnetar. 
The ejecta shell should decelerate in the long-lasting wind 
before it can hit the ion tail. A magnetized shock is likely to 
form in the wind instead of the ion tail, because the magneti- 
zation should be very low for the tail. Synchrotron maser can 
be produced as the blast wave sweeps the cold wind, which 
is illustrated in Figure 7. 
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Figure 7 The schematic picture of Beloborodov model. (Adapted from 
[219]. Reproduced by permission of the AAS.) 
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The evolution of I also has two stages separated by the 
crossing radius Rdec ~ 2T ?cót, 


ite) r < Raec, 


1/2 
2 (2c&E 
s ee 


where Lw, lI is the luminosity, Lorentz factor of the un- 
shocked wind and Lr = &/6t is the flare luminosity. Since 
the magnetization of the wind c, > 1, the peak frequency of 
maser emission is 


(77) 


r> Raec, 


~ 1/2 
To 3 e (=) r (78) 


yp = 3 f= DL. 
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Substituting r with t and using Eq. (77), the time evolution of 


Vpk is then 
tol tos, tobs < fo, 
Vpk = Vo x o obs T obs o (79) 
(to / tops) > tobs > to, 
where 
3/4 E 
= eL - Ly) (z) p 
Ime? E SPATS ^ ^ SNAG 54 (20 
T 1 ms 
to ~ — = —— ôt3. (80) 


Different from Metzger et al. (2019) scenario, He found that 
ICS is not important so that this intrinsic vpk is responsible for 
the observed peak frequency. Also, Eq. (79) predicts a differ- 
ent downward drifting rate compared to Eq. (76). These dis- 
tinctions originate from different physical conditions of up- 
stream medium. Note that due to high magnetization of the 
wind, a very low maser efficiency is also expected in this s- 
cenario, Sppp/& ~ 1077 rz? ~ 1075 — 1075. Moreover, FRB 
emission is almost 100% linearly polarized with polarization 
direction being aligned to the magnetar rotation axis since 
the helical wind is nearly toroidal. This is also disfavored by 
the observation of FRB 180301. A bubble structure can be 
formed behind the slow ion tail (Figure 7), and a strong op- 
tical flash of duration ~ 1s is predicted as the shock sweeps 
the bubble. However, no optical emission has been detected 
for FRB 200428 and all other events yet. 


6.3 A close-in model of Kumar & Bošnjak (2020) and Lu 
et al. (2020) 


Different from the previous two scenarios, the radiation 
mechanism in this model is coherent curvature radiation from 
bunches in the magnetosphere. A sudden disturbance (e.g., 
crustal quake) of the magnetar can launch Alfvén waves 
which moves along the field lines. In the polar cap region 
(i.e., close to the stellar surface), these Alfvén waves prop- 
agate outward where the density is decreasing. At a certain 
radius R,, the plasma density is too low to support the electric 
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current required by the wave. This means charge starvation 
occurs at R > Re. Strong parallel electric field forms and 
accelerates electron-positron pairs. The counter-streaming of 
pairs leads to two-stream instability and then bunch forma- 
tion. Afterwards, FRBs are produced by the coherent cur- 
vature radiation from bunches. This model is illustrated in 
Figure 8. 


Figure8 The schematic picture of Kumar & Bošnjak 2020 model. (Adapt- 
ed from [220]). 


Assuming the Alfvén waves have an amplitude óB and 
wavelength 4, the critical density as a function of R is giv- 
en by Ampere's equation, 
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where L represents the component perpendicular to the mag- 
netic field and R, is the radius of the magnetar. A paral- 
lel electric field develops in the charge starvation region and 
electrons and positrons are separated. According to the dis- 
cussion in Section 5.1.3, bunches of lengthscale l) ~ c/wp 
can form via two-stream instability. The transverse size of 
the coherent region is of order /, ~ R/y, therefore the num- 
ber of particles in the coherent region is 


Neon ~ nel (RI y^, (82) 


The total number of observed particles is larger by a factor 
(R/y*)/ lj, so the power of coherent curvature radiation is 
2e? cR yn? 


Po = Na RIY pe = s (83) 


where the single particle power p, = 2y*e?c/3p? is used. The 

isotropic luminosity of FRB is then 

16? cR^y?n?l, 
3p? ` 


Due to the uncertainties in these parameters, FRB luminosity 


Liso = 8y^P. = (84) 


could span several orders of magnitude. The magnetosphere- 
origin models have been applied to FRB 200428 and the 
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association with XRBs is also expected [124, 214, 225, 
226]. The derived spectra from CHIME and STARE2 
bursts can be explained [214]. Using the luminosity L ~ 
3 x 1038 erg s^!, the emission radius can be constrained as 
R./Rx ~ 20(6B10/A, 4)9/!! [124]. Note that these models pre- 
dict that X-rays should appear earlier than an FRB, since the 
FRB is not yet produced until Alfvén waves reach R.. How- 
ever, observationally there is a time lag of the X-ray peak 
respect to the radio peak, which is about several milliseconds 
[22, 23]. Further work needs to be done to explain this time 
lag as well as the non-thermal component of X-ray emission 
[22]. Also, the radio pulses should be able to break out from 
XRB fireballs [227]. 


6.4 A close-in model of Dai et al. (2016) and Dai (2020) 


Based on coherent curvature emission, Dai et al. (2016) pro- 
posed a different model, in which a repeating FRB source 
could arise from a strongly magnetized NS encountering an 
extragalactic asteroid belt (EAB) around a stellar-mass object 
[171]. The central source could either be a normal pulsar or 
a magnetar and the FRB emission region is also close to the 
stellar surface. In this model, an FRB could be produced by 
a collision of the NS with an asteroid of mass ~ 10" g to 
~ 10?! g in the EAB. Following Colgate & Petscheck [228], 
Dai et al. (2016) analyzed the collision physics [171]. It is 
assumed that an asteroid as a solid body falls freely in the 
gravitational field of the NS. This asteroid is originally ap- 
proximated as a sphere. It will first be distorted tidally by 
the NS at some break-up radius and subsequently elongated 
in the radial direction and compressed in the transverse direc- 
tion. The timescale of such a bar-shaped asteroid accreted on 
the NS surface is estimated by 


At = 1.6m;/? ms, (85) 


where and hereafter mig = m/10!8 g, the NS mass Mys = 
1.4Mo, and the asteroidal tensile strength and original mass 
density are taken for iron-nickel matter (see Equation 2 of 
[171]). The average rate of gravitational energy release near 
the stellar surface during At is approximated by 


Eg = GmMws (Rus At) = 1.2 x 10m? ergs |, — (86) 


where and hereafter the NS radius Rys = 106 cm is adopted. 
These simple estimates of At and Eg are well consistent with 
the observations of FRBs. This is why asteroid-NS collisions 
have been proposed as an origin model of FRBs [171,229]. 
Further, they also explored the radiation physics during an 
asteroid-neutron star collision in detail and found that an elec- 
tric field induced outside of the asteroid has such a strong 
component parallel to the stellar magnetic field that elec- 
trons are torn off the asteroidal surface and accelerated to 
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ultra-relativistic energies instantaneously [171]. Subsequent 
movement of these electrons along magnetic field lines will 
cause coherent curvature radiation. From Equation (15) of 
[171], the isotropic-equivalent emission luminosity is esti- 
mated by 


Liso ~ 2.6 X 10 (fo.e) mip pay ergs ^, — (87) 


where the beaming factor f is introduced based on the emis- 
sion beaming geometry [230], pc is the curvature radius of a 
magnetic field line near the stellar surface in units of 10° cm, 
u30 is the NS magnetic dipole moment in units of 10°° G cm?, 
and the other parameters are taken for an iron-nickel asteroid. 

This model was shown to be able to interpret many aspect- 
S of current observations under reasonable conditions. First, 
Dai et al. (2016) suggested that this model can not only ex- 
plain the observed emission frequency, duration and luminos- 
ity but also account for the repeating rate of FRB 121102 for 
the EAB's mass ~ 10 times the earth mass [171]. The down- 
ward drifting rate and linear polarization of FRB 121102 is 
well consistent with model prediction [231 ]. 

Second, the periodical activity of the repeating FRB 
180916 can be understood. Repeating bursts from this source 
show a ~ 16 day period with an active phase of ~ 4.0 days 
and a broken power-law form in differential energy distribu- 
tion. Dai & Zhong (2020) proposed that FRB 180916-like pe- 
riodical FRBs could provide a unique probe of EABs [232], 
following the model in which repeating FRBs arise from an 
old-aged, slowly spinning, strongly magnetized NS traveling 
through an EAB around another stellar-mass object [171]. 
These two objects form a binary and thus the observed period 
is in fact the orbital period. Furthermore, the EAB's physical 
properties were constrained and it was found that (1) the out- 
er radius of the EAB is at least an order of magnitude smaller 
than that of its analog in the solar system, (2) the differential 
size distribution of the EAB's asteroids at small diameters 
(large diameters) is shallower (steeper) than that of solar sys- 
tem small objects, and (3) the two belts have a comparable 
mass. 

Third, Dai (2020) proposed a specific model for FRB 
200428 and its associated XRB, in which a magnetar en- 
counters an asteroid of mass ~ 107° g [230]. This asteroid 
in the magnetar's gravitational field is first possibly disrupt- 
ed tidally into a great number of fragments at radius ~ a few 
times 10! cm, and then slowed down around the magnetic 
interaction radius by an ultra-strong magnetic field and in the 
meantime two major fragments of mass ~ 10!" g that cross 
magnetic field lines give rise to two pulses of FRB 200428, 
as shown in [171]. Dai (2020) also argued that the whole 
asteroid is eventually accreted onto the poles along magnet- 
ic field lines, colliding with the stellar surface, generating a 
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photon-e* pair fireball trapped initially in the stellar magne- 
tosphere, and further leading to an XRB [230]. This model 
was shown to explain all of the observed features of FRB 
200428 and its associated XRB self-consistently. In addition, 
this model can interpret the time delay ~ 3 — 6ms of the 
XRB’s two peaks with respect to FRB 200428's two puls- 
es, which was indicated by the observations of the Insight- 
HXMT and INTEGRAL satellites. This time delay disfavors 
a few models, e.g., a starquake-induced close-in model [124], 
a synchrotron-maser far-away model [215], and also the other 
magnetar-asteroid model of Geng et al. (2020) [233]. Finally, 
this model is well consistent with the constraints on bursting 
sites inferred from both the Insight-HXMT’s detection of the 
time interval between the XRB's two peaks and the time in- 
terval between FRB 200428's two pulses. 


7 FRBs as Cosmological Probes 


In the near future, an FRB sample enlarged by orders of mag- 
nitude will be obtained, with a considerable fraction being lo- 
calized and covering a wide range of redshift. It can be used 
to probe cosmological parameters such as the energy density 
of matter, dark energy, cosmic curvature, the baryon fraction 
of the IGM, the dark energy equation of state and the reion- 
ization histories for helium and hydrogen [234]. This can be 
realized by means of the DM(z) relation given by Eq.(9). 


7.1 Cosmic Baryons 


For the current universe, the significant discrepancy between 
the detected baryonic matter and the theoretically predict- 
ed baryonic matter is called the missing baryon problem 
[235, 236] . It is difficult for us to detect baryonic matter 
in the diffuse gas of the universe. Big Bang nucleosynthe- 
sis and Cosmic Microwave Background (CMB) can be used 
to constrain the density of cosmic baryons [237, 238]. There 
are many works that use simulated FRB data samples to con- 
strain fiom [32,239-241] and try to alleviate the missing bary- 
on problem [242-244]. Due to the large DMs of the cosmo- 
logical FRBs, the localized events can be used to measure 
the baryonic matter of the universe. Once the DM value and 
the redshift of the localized FRBs are obtained, the DMjgm 
value contributed by the intergalactic medium can be calcu- 
lated. Only DMjgm is related to baryon density, the value of 
which can be determined by combining redshift. Macquart 
et al. (2020) used the five ASKAP localized FRBs to derive 
a cosmic baryon density of Q, = 00575. 
shows a great consistency with the Big Bang nucleosynthesis 
and CMB measurements, indicating that it is feasible to mea- 
sure the density of cosmic baryons using localized FRBs (see 


and their result 
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Figure 3 of [82]). 


7.2 Dark Energy 


In cosmology, dark energy is devoted to explain the accelerat- 
ing expansion of the universe. The extragalactic FRBs can be 
used as a cosmic probe to study the evolution of the universe 
and constrain cosmological parameters. A Combination of 
FRB with other cosmological probes will be a powerful tool 
in constraining the dark energy equation of state [245-250]. 
The Friedmann equation can be written as [251]: 


H(z) 
HG 


= Qmu(1 +z)? + Qpe exp IET  oXz)dIn(1 +z) 
+Q + zy), (88) 


where Oy and Opg are the matter and dark energy density of 
the universe, € is the spatial curvature density. w(z) is the 
dark energy equation of state. The constraints on dark energy 
can be obtained using the DM(z) relation. 

Zhou et al. (2014) constrained the equation of state 
of dark energy using type Ia supernovae (SNe Ia), simu- 
lated FRBs and BAO data in the optimistic case that in 
each narrow-redshift-bin tens of FRBs have been detected 
[245]. FRB can also be associated with other cosmological 
probes to constrain the dark energy equation of state, such as 
GRB/FRB system [247], lensed FRBs/CMB/SNe Ia system 
[246], GW/FRB system [248], CMB+GW-+EFRB [250] and 
FRB/standard distance probes [249]. However, the effects of 
DMs contributed by host galaxies and fluctuations of IGM 
should be treated properly. 


7.3 Proper Distance 


The proper distance dp is the length from the source to the 
observer along the light of sight at a specific moment of cos- 
mological time. In the frame of the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric, it is expressed as 


r d / 
dy(r)-a i E E (89) 
0 N1- Kr? 


where a is the scale factor, r is the co-moving coordinate of 
the source and K is the curvature parameter. 

The proper distance is difficult to measure because it de- 
pends on the cosmic curvature. FRB is a good probe that 
is independent of the cosmic curvature and record infor- 
mation about the expansion of the universe. The proper 
distance-redshift relation can be derived from the DMs of 
FRBs with measured redshifts. Yu & Wang (2017) proposed 
that dy — z relation can be derived from 500 mocked FRBs 
with (DMiay, z) data through Monte Carlo simulation [252]. 
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Figure 9 shows the ability of FRBs in measuring proper dis- 
tance from simulations. 500 FRBs with DM and redshift 
measurements can tightly constrain cosmic proper distance. 
This dy —z relation can constrain the cosmic curvature as well. 
Of course, this depends on the number of FRB detections in 
the near future. With the advancement of observational tech- 
nique, this will become an effective method to measure the 
inherent distance. 
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Figure 9 Top panel shows the binned mock DMjgm data with lo error 
bars. Bottom panel shows the dp(z) function with its 1o confidence region 
derived from Gaussian Process method and its theoretical function which is 
shown by red line. Qk = 0 is assumed. (See [252] for more details.) 


7.4 Hubble Parameter H(z) 


FRBs with redshift measurements can be used to measure 
Hubble parameter H(z), which is a very important param- 
eter in cosmology. In flat ACDM cosmology, H(z) can be 
expresses as: 


A(z) = Ho VOR + Oyw(1 t zy, (90) 


where Q; is the vacuum energy density fraction. 

Wu et al. (2020) proposed a new model-independent 
method to calculate the Hubble parameter H(z) [253]. The 
averaged redshifts (z) and dispersion measures (DMjgm) of 
FRBs are known. Differentiating equation (9), the Hubble 
parameter can be expressed as 


HO = ORF) (91) 


&nGm, ADMicm’ 


where Az and ADMjcp are the differences of z and DMjgm 
between two adjacent bins, respectively. Figure 10 shows the 
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simulation result. The error of derived H(z) is about 696. The 
DMs and redshifts of 500 simulated FRBs are used to mea- 
sure H(z), and the derived results are consistent with the theo- 
retical values of H(z). As more localized FRBs are observed, 
the calculation of Hubble parameters that rely on redshift will 
become possible and provide useful probes for cosmology. 
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Figure 10 A(z) and 1c errors derived from 500 mock FRBs are shown as 
blue dots. The red line represents the theoretical H(z) function. The devia- 
tion between the simulated value of H(z) and the theoretical one is about 696 
at z = 2.4. (Adapted from [253].) 


7.5 Dark Matter 


Dark matter is abundant in the universe. Because it does not 
interact with electromagnetic fields, it is difficult to detect it 
directly. Massive compact halo objects (MACHOs) are the 
possible candidate for dark matter, and the fraction of dark 
matter fom in MACHO has received a lot of attention. Strong 
lensing of FRB by MACHO can be used to constrain fpm. 
When there is MACHO in the LOS, due to gravitational lens- 
ing effect, two images and the corresponding time interval 
will appear. For a MACHO with mass M, as a point lens, the 
Einstein radius is 


GM, Drs 
c? Ds Dr’ 


On = 2 (92) 
where Ds, Dz and Dzs are the distances to the source, to the 
lens and between the source and the lens, respectively. The 
time delay between the two separated images is [254,255] 


N CN oe | 
2 Vy+4-y 
where y = 8/0g and £ is the impact parameter. 

Muñoz et al. (2016) proposed that the lensed FRBs with 
double peaks can be used to constrain the fraction of dark 


matter in MACHOs [254]. They constrained fp, in MA- 
CHOs to fpm < 0.08 for Mi, = 20Mo if there are no FRBs mi- 


|. 4GM, 
= 


At (1 + zz) (93) 
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crolensed. Wang & Wang (2018) proposed to probe the com- 
pact dark matter with FRBs gravitational lensed by MACHO 
binaries [256]. They concluded that fpm in MACHO is con- 
strained to « 0.001 when there were no search of multiple- 
peaked FRBs for time intervals larger than 1 ms and flux ratio 
less than 10?. These methods of using FRBs to find dark mat- 
ter in MACHO is widely studied later [257-261 ]. 


7.6 Hubble Constant 


The large event rate, small temporal duration and cosmologi- 
cal origin of FRBs suggests that FRBs can be a clean cosmo- 
logical probe. Li et al. (2018) proposed that strongly lensed 
repeating FRBs can constrain Hubble constant Ho [262]. The 
accurate time intervals between different images of lensed 
FRBs can establish connections with Hubble constant, which 
is a precision probe of the universe. The Hubble constant Ho 
can be constrained with ~ 0.91% uncertainty by simulated 
10 lensed FRB systems within the flat ACDM model. Com- 
pared with other possible constrains on Ho, lensed FRB has 
obvious advantages in ascendancy. In addition, cosmologi- 
cal curvature can also be derived using lensed FRB, which is 
independent of model [262]. The gravitational lensed FRB- 
s can be used in various cosmological research in the future 
[263]. 


7.7 Weak Equivalence Principle 


Einstein's Equivalence Principle (EEP) is the equivalence of 
gravitational and inertial mass, which is the foundation of the 
General Relativity. The verification of the EEP comes from 
the parameterized post-Newtonian (PPN) parameters, such as 
the parameter y. Calculating the difference of the y value be- 
tween different particles ^y can constrain the accuracy of the 
EEP. Now, the extragalactic transients have been used to the 
constrains on EEP or Einstein's Weak Equivalence Principle 
(WEP) [264-268]. 

Wei et al. (2015) proposed that FRBs with cosmologi- 
cal origin can be a good candidate for constraining the EEP 
through the association with GRB [269]. They obtained that 
the upper limit of the difference of the PPN parameter y is 
[y(1.23GHz) - y(1.45GHz)] < 4.36 x 107°, which can be 
expressed as 


-1 

Ay = 71-72 < (Atos = Atom) (29) ta (5): (94) 
where Afops is the observed time delay between two differ- 
ent energy bands, Atpy is the time delay contribution from 
the dispersion, Myw is the Milky Way mass, d is the lumi- 
nosity distance from the source to the earth and b is the im- 
pact parameter of light rays relative to the Milky Way center. 
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This is a strong restriction on the EEP, which is more accu- 
rate than the restriction using SN and GRB. Using the above 
method, Tingay & Kalpan (2016) used the first localized FR- 
B 150418 alone to limit the EEP [270]. They obtained a limit 
of Ay < 1-2 x 107°, which is an order of magnitude better 
than the previous constraint [269]. 


7.8 Contamination from DMpos and IGM inhomogene- 


ity 
7.8.1 DMhost 


In order to know the exact value of DMjgm, the other 
contributions to the observed DM should be well under- 
stood. In Eq. (7), DMmw represents the total contribution 
from the Milky Way, which consists of two components, 
DMyw = DMgaisk + DMhato. The disk component DM, can 
be subtracted using the NE2001 model [33] or YMW16 mod- 
el [34]. DMhalo is the part contributed by the halo of Milky 
Way, of which the typical value is 30 — 80 pc cm^? [85,271]. 
The source contribution DMgource depends on the environ- 
ment of the progenitor. If FRBs are produced by the rem- 
nants from collapses of massive stars, the values of DMsource 
can be large [272]. Differently, if FRBs are produced by BNS 
mergers [273], DMgource is usually small [98, 274]. Once the 
central engine is known, DMgource can be derived analytically 
[98, 272, 275]. 


The contribution by host galaxies DMhost is poorly known. 
Xu & Han (2015) estimated the DM contributions by differ- 
ent types of galaxies using the scaled model of NE2001 [276]. 
Methods of constraining DMnhost statistically have been dis- 
cussed [277,278]. Moreover, the IlustrisTNG simulation is 
used to derive the distributions of DMpost [279]. Figure 11 
gives DMnhost at different redshifts. For repeating FRBs, the 
median DMpost are 35(1 + z)! 5 and 96(1 + z)°*? pe cm^? for 
FRBs like FRB 121102 and FRB 180916, respectively. The 
median of DMhost is about 30 — 70 pc cm ^? for non-repeating 
FRBs in the redshift range z = 0.1 — 1.5. In this case, the 
evolution of the median DM, can be fitted as 33(1 + z)°*4 
pe cm^?. The distributions of DMpost of repeating and non- 
repeating FRBs can be well fitted with the log-normal func- 
tion 
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Figure 11 The distributions of DMhost at different redshifts. The red lines 
are the best-fitting results using log-normal distribution. The blue shaded 
region in z = 0.2 panel is the DMpost for FRB 121102. (See [279] for more 
details.) 


7.8.2 IGM inhomogeneity 


McQuinn (2014) considered the effect of inhomogeneity with 
three models for halo gas profile of the ionized baryons [280]. 
Dolag et al. (2015) and Pol et al. (2019) studied DMjgm with 
different cosmological simulations at low-redshift (z < 2) uni- 
verse [271,281,282]. Efforts of tacking IGM inhomogeneity 
using the FRB data have been made [283, 284]. Jaroszyns- 
ki (2019) used the Illustris simulation to estimate the DMjGw 
and its scatters in the z « 5 universe [285]. Zhang et al. (2020) 
and Takahashi et al. (2020) used the IllustrisTNG simulation 
to realistically estimate the DMjgm up to z ~ 9 [279,286]. At 
low redshifts, the distributions of DMjgm are shown in Figure 
12. They are well fitted by [82] 


A^? — Co)? 
picm(A) = AA? ge N A»0. (96) 
20° Os 


where A = DMigm/{(DMjcm), and f is related to the inner 
density profile of gas (p œ r ^) in halos. opm is an effective 


Xiao, Wang & Dai, et al. Sci. China-Phys. Mech. Astron. 


standard deviation. C, can affect the horizontal position. 
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Figure 12 The distributions of DMjgm™ at different redshifts. The red lines 
are the best-fitting results using Eq. (96). (Adapted from [279].) 


In order to obtain the reliable cosmological constraint from 
FRBs, the distributions of DMpost and DMjay must be prop- 
erly handled [82]. The potential of FRBs as cosmological 
probes is huge, however, it is currently subject to limited sam- 
ple of localized bursts. With more events coming soon, FRB 
could become a very powerful probe. 


7.9 Constraints on photon mass 


The rest mass of a photon is supposed to be strictly zero and 
the velocity of the photon is the speed of light. It is very im- 
portant to test the correctness of this hypothesis experimen- 
tally, because it involves many fundamental physical prob- 
lems. The arrival time of the lower frequency pulse of an 
FRB is later than the higher one due to plasma effect. The 
frequency dependence of the pulse arrival time of the FR- 
B follows the v? law, which can constrain the photon mass 
[287-289]. Considering two photons with different frequen- 
cies, the high frequency vp and the low frequency v. The 
photon mass m, can be constrained as 


2HoAL, ve 
- , (97) 


-2 
my, = hc | —— 
lc - d H(z) 
where 
(1  z)?dz 


Hos | (98) 
Qa z) + Qn 


and At,,, is the arrival time interval of the low frequency and 
the high frequency photons. 

Based on the observations of FRB 150418 at z — 0.492, 
Wu et al. (2016) found an upper limit of m, < 5.2 x 
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107? g, which is a strict constraint on the photon mass 
[287]. The plasma effect and scattering due to multi-path 
propagation through the intervening material should be con- 
sidered. Bonetti et al. (2016) derived the photon mass 
my, € 3.2x 107?? ke by considering the frequency depen- 
dence of dispersion [288]. Shao & Zhang (2017) developed a 
Bayesian framework to constrain the photon mass with a cat- 
alog of FRBs [289]. With the growth of FRB population, the 
constraint on the photon mass will be promising using this 
method. 


7.10 Reionization history 


The transition from neutral to ionized universe is called reion- 
ization. Usually, it is believed that the reionization of hydro- 
gen occurs between 6 « z « 12. For helium the range is 2 « 
z « 6 [290]. When reionization happens is still unresolved 
in the modern cosmology. Extragalactic FRBs can be a good 
probe for tracing the history of the H and He II reionization 
using high DM contributed by IGM [32, 150,242,291]. Ac- 
cording to Eq.(9), the ionization fractions of intergalactic hy- 
drogen X. and helium Xe He can be derived when the red- 
shifts and DMjay of FRBs are measured. 

Zheng et al. (2014) proposed that FRBs can be the po- 
tential probe to study the history of He II reionization [242]. 
During the epoch of the He II reionization at redshift z ~ 3, 
the derivation of DM value decreased 8%. This indicates that 
DM will changes across the He II reionization epoch. It is 
clearly that a large FRB population is needed to study the He 
II reionization history. Caleb et al. (2019) simulated the FR- 
B samples, and tracked the evolution history of reionization 
using the simulations of DM values of FRBs. They obtained 
the number of FRBs needed to distinguish reionization oc- 
curred at z = 3 or z = 6 through simulations. They concluded 
that more than 1100 FRBs were needed when the intergalac- 
tic medium IGM is homogeneous and more than 1600 FRBs 
needed when the IGM is inhomogeneous [150]. The joint ob- 
servation of the next generation radio telescope will observe 
enough FRBs to probe the history of the He II reionization. 


7.11 Intergalactic turbulence 


The intergalactic turbulence is very common in the universe, 
which is related to the formation of large-scale structure of 
the universe. The observations and statistical analysis of the 
intergalactic turbulence are very useful in astronomical re- 
search [292]. The high DM of FRBs is mainly contributed by 
IGM, so the dispersion disturbance of DMjgm can reflect the 
density fluctuation in the IGM. Thus, an FRB can be used to 
probe intergalactic turbulence. Recently, Xu & Zhang (2020) 
first proposed a statistical measurement of the intergalactic 
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turbulence by using a population of FRBs [292]. The struc- 
ture function of DMs can represent the multi-scale electron 
number density fluctuations from the intergalactic turbulence. 
Xu & Zhang (2020) calculated the structure formation value 
using the published FRB data, which has large amplitude and 
a Kolmogorov power-law scaling. They concluded that the 
intergalactic turbulence extended outward to 100 Mpc with a 
Kolmogorov power spectrum. The extensive application of 
using fast radio bursts to study the intergalactic turbulence 
depends on the accumulation of more FRB data in the future. 


8 Summary and Future Prospects 


The FRB field is very energetic and rapidly developing. Im- 
portant findings keep coming quickly. The current status and 
major developments in understanding their inherent physics 
are summarized here. 


() It is controversial whether all FRBs repeat. Sever- 
al statistical properties of repeaters and apparent non- 
repeaters seem different. However, no clear sign of di- 
chotomy has been found. It has been claimed that these 
differences could be owing to the diversity in repeating 
modes. The other viewpoint that some authors prefer is 
that nature is diverse, and two populations could have 
different physical origins. The existence of genuine- 
ly non-repeating FRBs could be finally confirmed if an 
FRB from a cataclysmic event (e.g., BNS merger) is ob- 
served. The large offset and host galaxy property of FR- 
B 180924 and FRB 190523 seem to favor this hypothe- 
sis. As may be imagined, it will be a big breakthrough 
if FRB is associated with gravitational-wave signal and 
other EM counterparts. 


(ii) The radiation mechanism of FRBs remains unknown. 
We have introduced five mechanisms, which do not ex- 
plicitly define the whole mechanism. Currently, the co- 
herent curvature emission and synchrotron maser emis- 
sion from magnetized shocks are in the leading position. 
These two have been studied in great detail and applied 
in the case of FRB 121102 and FRB 200428. Under a 
few assumptions, both of them can explain the observa- 
tional data well. Their predictions are different and ver- 
ifiable in the future. Some pieces of emerging evidence, 
such as the FAST detection of diverse polarization an- 
gle swings from FRB 180301 [69] support the former 
emission mechanism, i.e., a pulsar's magnetospheric o- 
rigin. 


(iii) The magnetar origin of FRBs has been established 
through FRB 200428. However, the emission site has 
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not yet been determined. Both close-in and far-away 
models can explain the observations well [216, 230]. 
(2020) reported a repetition of FRB 
200428 with double peaks again after long-time mon- 
itoring [21], which makes the situation more interest- 
ing. However, conclusive evidence has not yet been 
found, which might be polarization property. General- 
ly, far-away maser models predict that FRB is linearly 
polarized with a flat PA curve, while close-in models 
could have diverse PA swings according to the geome- 
try of magnetic field and observing direction. This will 
come to a conclusion once the polarization information 
of FRB 200428 is obtained from the subsequent repeat- 
ing bursts. 
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(iv) The size of the total FRB sample is still small. We 
definitely need more FRBs to precisely analyze the sta- 
tistical properties. Especially, the number of localized 
FRBs is very limited. Hundreds of localized FRBs are 
needed to do fine cosmology. In this sense, the popu- 
lation study and cosmological application of FRBs are 
still at an early stage, and plenty of work needs to be 
done in the future. With the upgrade of current instru- 
ments and the upcoming new facilities, the total number 
should reach the order of a thousand in a short time. 


Of course, certain open key questions in the FRB field re- 
main to be addressed over the next few years. Do genuinely 
non-repeating FRBs exist? What is the real ratio of repeaters? 
Do all repeaters have a periodical activity? What is the form 
of the intrinsic luminosity function and redshift distribution? 
Is there any other physical-motivated criterion for classify- 
ing FRB? Except for magnetars, could any other progenitors 
produce FRBs? Are there multi-wavelength counterparts or 
multi-messenger signals accompanying FRBs? These ques- 
tions are expected to be answered sooner or later, with future 
observations performed by high-sensitivity radio telescopes, 
such as FAST and SKA. 
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